// RNG works on C++ compilers (e.g., Microsoft Visual C++) w 64-byte integer
inline double RNG ( _int64 *pSeed)
{
  *pSeed = (*pSeed * 16807) % 2147483647;
  return (double) *pSeed / 2147483647;
}
